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1. INTRODUCTION 

Wheel mobile robot (WMR) are simple in structure, certainty, and high energy efficiency [1]. The 
WMR have many applications in various fields such as transportation, military, medical, and rescue. 
Underactuated WMR which is popular among mobile robots, is a non-linear subject with multi-input, 
multiple-output, and has nonholonomic constraints. As a result, this WMR has attracted the interest of many 
researchers working on modeling and control algorithms [1]-[3]. 

In recent decades, many literatures have been published with a variety of control methods for tracking 
control of WMR. Among these, the sliding mode control (SMC) [4]-[10] is one of the most popular technique 
which is employed for WMR to cancel or reduce the effect of unknown components to the working of the 
system. The control schemes in [4], [5] are traditional SMC for closed-loop dynamics. The cascade control 
structure provides a good control quality and becomes a favorable structure for several current researches. 
Also, cascade control structure is introduced in [6] for tracking control of WMR. The kinematic controller is 
designed based on Lyapunov function and the dynamic control loop is stabilized by SMC. The presented 
control algorithm guarantees the fast speed tracking and eliminates the kinematic uncertainties. According to 
Goswami and Padhy [7], the SMC is employed for both inner-loop and outer-loop to cancel the effect of 
disturbances. In this work, the wheel slips are integrated into the system model and this becomes popular in 
later works. According to Zhai and Song [8], the nonsingular terminal SMC without any constraints is 
developed for WMR with external disturbances and inertia uncertainties considering. The proposed control 
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scheme ensures the finite-time convergence but the tracking error only converges to a small region. Although 
all of the above SMCs provide good performances for the tracking control problems of the WMR, a common 
drawback of SMC, chattering phenomena, still exists. For this reason, many researchers are constantly studying 
and proposing algorithms to improve the performances of the SMC systems. 

Parallel with the SMC, adaptive control techniques [11]—[25] are also employed to solve the problems 
of the WMR systems. According to Amer et al. [13], the backstepping control is used for kinematic loop and 
fuzzy SMC for dynamic loop to achieve the robustness and the stability. The fuzzy logic system (FLS) is 
utilized to tune the adaptation gain which increases the adaptability and chattering free. However, despite the 
reasonability of the theory, the demonstrated result in this work is unclear. Another approach of the adaptive 
control scheme which employed the neural network (NN) controller based on the reinforcement learning 
method is proposed in [20]. In this work, NNs are trained to address the uncertainties, the slipping and the 
skidding components which affect the tracking quality of the system. With this structure, the system 
performances are improved countably; however, the complexity is one of the disadvantages of this scheme. 
Designing an adaptive control in combination with a disturbance observer is also one of the most effective 
ways to compensate the effect of disturbance [22], [23]. Especially, according to Li et al. [22], the cascade 
structure which incorporates a nonlinear disturbance observer with the extended Kalman filter is used to avoid 
the influence of disturbance. Although the wheel slipping can be compensated, the overshoot of the 
performance is still too high. Meanwhile, the strategy in [23] is also the two-loop structure with a disturbance 
observer in the second loop. With the addition of the observer, the feedforward compensation for the lumped 
disturbance is obtained, but the settling time in the result is large. Besides, to improve the control quality, [24] 
also proposes a control scheme which combines the sliding-mode control and disturbance observer. In this 
case, two observers are executed to eliminate the effect of the disturbance. However, this structure increases 
the number of complex calculations tremendously. 

This paper introduces an adaptive SMC for the purpose of tracking control of a WMR in the presence 
of uncertainties, input disturbances, and wheel slips. The establishment of the state space model of the system 
is primarily based on the kinematic and dynamic models. Subsequently, the state space model is stabilized 
through the implementation of the adaptive SMC scheme. The unidentified elements within the system are 
approximated through the utilization of an observer, which subsequently transmits the information to the 
controller in order to adjust for the impact of any disruptions. The observer's convergence is guaranteed through 
the utilization of Hurwitz criteria, while the stability of the entire system is established on the basis of the 
Lyapunov theory. The advantages of the proposed scheme are stated as: i) the control system is designed with 
only one loop instead of two loops as usual. This structure reduces the complexity in both study and 
implementation and ii) the control signal is free from chattering phenomena of the SMC which gives the 
reduction of the cost of energy, extending the acuator’s lifespan. Finally, the comparative simulations are done 
to demonstrate the advantages of the proposed adaptive SMC in comparison with the existing work. 


2. MATHEMATICAL MODEL 

In the article, the study focuses on the WMR concerning the robot's motion on a horizontal surface. 
The body of the system incorporates two rear-wheel drive wheels for propulsion and a front steering wheel. 
The movement frames of the WMR are illustrated in Figure 1. 


YM 


M 


Figure 1. Wheel mobile robot model 


Let Yr, Y, are the longitudinal slip factor of the left and right driving wheels, n is the lateral slip factor 
along the wheel shaft. Angular velocity of the right and left wheel are g, È, respectively. Linear forward 
velocity in perpendicular direction of the wheel shaft and angular velocity are calculated as (1) [12]: 
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9 = Ortha) + rth, ai rnd) + Lae (1) 


The mathematical model of the system is built based on the kinematic and dynamic equations. The kinematic 
and dynamic equations of the WMR are described as [12]: 


žm = V cos 0 — ù sin 0 


Yu = V sinð +)cos08 (2) 
ĝ =w 
Mù + Bu+Qyt+Cn+Gi+ty=T (3) 
where: 
. . Mı M42 
v=[de gi)", yv=lYr Y], M = bn al M22 = M11, M21 = M12, 
2 2 2 2 22 2 
my, = me (+ 2) + (Ug + 2lp) + 2myr? + ly, My = mg (S-E) - 7 Ug + 2lp) 


= i 0 1 _ [21 Q = r a? r ~ r fl 
B = mz” Es al Q= 0, A! Qi2 = MG “(4 ar z) + zz (le + 2ip), C = m70 hi 
game [3] 
€ 2b l-1 
In which mg is WMR body mass; mw is each driving wheel mass; Jc is the inertial moment of the platform 
about the vertical axis through point G; Jp is the inertial moment of each wheel about its diameter axis and; Iw 
is the inertial moment of each wheel about its rotational axis. Formula (3) can be rewritten as (4) and (5): 
v=—-M"Bv+M'4*t-—M1(Q¥ + CH) + GÑ + Ta (4) 
ù = —-M7tBv + Mtt +d, (5) 
Where: d) = —MTI (QY + CH) + GÑ + Ta 


For tracking control, define target point as Z (xz, yz). The position error between point M (xm, Ym) and Z in 
OXY is as (6) [12]: 


ep = | = l cosð® sin j — Xm 


~ Lep2 —sin@ cos@ll¥z— Ym (6) 


In (6) expresses the tracking error of the control system. The derivative of (6) along with time under the 
condition of wheels slip and external disturbances is obtained as (7): 


ie Nee en e 
é, = a =hv +d, (7) 


where: 


ep2 r ep2 r ÝR-ÝL _ VRTYL ; 
_ (2 aa) =e +1); w ( 2b )ep2 2 cos@ sin 0) [xz 
h aa , dı res . . + S A 
Spar pir — (4) y —sin@ cos 8l yz 
ep1 n 
2b 2b 2b 


Define variables x, = €p; X2 = X, + Ax,, where A is a positive scalar. The time derivative of x, is: 


X2 = ži + Axı (8) 


To be convenient to the system model build, (8) be rewritten to cancel the high order derivative of x. 
Substituting (5) and (7) into (8) leads to: 


X = hý + İv + d; + A(hv + dy) 
=hM-1Bv +hM1t + hd, + hv +d, + Ahv + cd, (9) 
= Pv + Et + d3 
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in which: 
P,=—-hM1B, E=hM"1, d3=hd,+hv+d,+Ahv+dd, 
From (9) can be reformulated in the form of state variable. By substituting v from (7) into (9), x2 is written as (10): 
ž2 = Px —APx,+Et+d (10) 
where: 
P=P,h*=M"B, d=d; -hd 


From above calculation, the control model of the WMR will be established. The WMR state model can be 
expressed as (11): 


i = X — Àxı (11) 
žə = Px, —APx,+Et+d 


In (11), A is the positive constant with unknown exact value. For the good control performance, in 
this work, A is viewed as unknown component which needs to be estimated. In this sense, define disturbances 
Dı = Ax, and D, = d — APx,. The state model (11) becomes: 


fa = x, + Dı 


ty = Px, + Et + D, 2) 


3. ADAPTIVE SLIDING MODE TRACKING CONTROLLER DESIGN FOR WMR 

In this section, an adaptive SMC will be introduced for tracking control of the WMR. The control 
scheme is designed through two steps: the first one, the chattering free SMC is built based on the real 
disturbances. Next, the unknown disturbances will be estimated by a disturbance observer then the real 
disturbances will be replaced by estimated ones. The detail of each step is illustrated in the following 
subsections. 


3.1. Disturbance observer based stable sliding surface design 
Before designing the SMC, the sliding surface should be designed first. The sliding surface is defined 
as (13): 


S = ž + kox. +k x, (13) 


As the sliding surface is reached, s = $ = 0. Under this condition, replace (13) into (12) leads to the as (14): 


ma ae (14) 


Xq = —k2x3 — kyx, 


In (14), if the mismatched disturbance D;=0 then system states converge to the stable points. However, D1 is 
different from zero so from (14) the following result is obtained: 


ži + k2%1 + kx = Dı + k2D, (15) 


This means that the stability of x; can not be guaranteed because of the effect of the mismatched 
disturbance D;. To overcome this drawback, the nonlinear disturbance observer [25] will be used to estimate 
the mismatched disturbances. 


B = §,—a,(f, — x1) + x2 


; (16) 
> = —a, tanh(b(, — x1)) 


Where é; and & are estimated values of x; and D;, respectively; aj, a2, and b are positive scalars. 
Define new variables e} = €, — x1, €a = Dy — &4; eg = —eq — a, e,. Considering (14) and (16), the 
derivative of e; and e2 is as (17): 


ie = (17) 


é, = —a, tanh(be,) — D, — a,e, 
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The convergence of e; and ez as well as ea is given detail in [25]. Substituting (16) and (17) in to (15) 
and do some manipulations leading to as (18): 


Xt kaž + kyxy = bg + keg (18) 


Because the estimated error eg is bounded, it can be realized from (18) that if k; and k2 are Hurwitz then the 
system states are stable. In considering the disturbances, the new sliding surface is established as (19): 


S = ty + ko (x2 + &) + kıxı + & (19) 


3.2. Adaptive sliding mode control design 
After ensuring that the sliding surface (19) is stable, in this section, the reaching controller will be 
designed for closed loop stable purpose. Introduce the control law [25]: 


T=Q1(t, +72) (20) 
where: 

Ty = —PX_ — ky (xq + &) — kx, — & (21) 
and 

Tz tat, = —(e +k3 + u) sgn(s) (22) 


in which a > 0, € > |D, |, k; > a|D,|, and py is a positive scalar. Substituting (11) into (19): 
s = Px, + Et + Dy + ko (x2 + §2) + ky + & (23) 


The sliding surface (23) can be shortened by replacing the control law (20) with using (21) and (22) into (23). 
By this way, the sliding surface is obtained as (24): 


S=T2 + D, (24) 
The stability of the reaching controller is proven via Lyapunov theory. Choosing the Lyapunov function: 


V =-s's (25) 


2 


Next, the time derivative of V will be proven to be negative. The time derivative of V along with (21) and (22) 
is as (26): 


V =s's =s"(t, + D2) = sT (i, + aT, — at, + D2) 
= s’[-(e + k; + u) sgn(s) — at, + D3] 
= (—e|s| + s™D,) — (k3s7 + as™T2) — pIs| 
< —u|s| <0 (26) 


This implies that the introduced controller (20) guarantees the reaching to the sliding surface of the system 
states as time goes to finite. 


4. SIMULATION RESULT 
In order to verify the effectiveness of the proposed control scheme, the simulations are setup in 
MATLAB/Simulink for a WMR with the following parameters: 


Mg = 10kg, m, = 2kg, ly = 0.1kgm?, In = 0.05kgm?, Iç = 4kgm?, b = 0.3m, r = 0.15m, 
a=0.2m 


Parameters of the controller are given by: à = 10,k, = 10, k, = 10,2 = 100,u = 10,£ =5. 
The simulations are executed in considering the input disturbances as (27): 


Tq = [2 + sin(0.1t) 1+cos(0.1t)]" (27) 
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and the longitudinal and lateral slip factor can be chosen as (28): 


[0, 0,0]? (m/s) t < 2(s) 


[0.3sin(2t), 0.3cos(2t), 0.2]? (m/s)t = 2(s) (28) 


eV)’ = { 
The WMR is controlled to move with two kinds of reference trajectories: circular shape and trifolium 
shape. To demonstrate the effectiveness of the proposed adaptive sliding mode control (ASMC), this paper 
simulates the design algorithm along with the represented method in [6]. The simulation results for each case 
are illustrated in the following. 
- Case 1: circle shape with the following equations of reference trajectory: 


Xz =5cos (= t) „Yz = 5sin (= t) (29) 


The initial conditions of the mobile robot are given by: [x,(0), yy (0), 0 (0)] = [5; —3.2; 2/2]. 

The trajectory tracking of the ASMC and Department of Statistics Malaysia (DOSM) is depicted in 
Figure 2. It can be seen from Figure 2 that the responses of both ASMC and DOSM are good with similar 
trajectory at steady state time. However, at initial time, the response of the DOSM takes a long time to track 
the reference curve and has little oscillation at 3 s where the disturbances act on the system. This result is 
expressed thoroughly in Figure 3(a) which presents the tracking error in x-axis and y-axis. In Figure 3(a), the 
tracking errors of two controllers are almost the same in the y-axis with nearly zero steady state error. In the 
x-axis, the maximum tracking error of the proposed ASMC is about 0.05 m while this is bigger for DOSM 
(about 0.65 m). The advantages of the proposed ASMC and the DOSM are illustrated clearly in Figure 3(b), 
the control torques. In Figure 3(b), the control torques of the ASMC have a big pulse at initial, but they are 
smooth during the working time. Meanwhile, the control torques generated by the DOSM scheme are affected 
by chattering phenomena of the SMC, so they are highly oscillated in both transient time and steady state. This 
chattering is not good for power electronic devices as well as the mechanical systems. Also, at the time acting 
on the system, the control torques of the DOSM are too higher than those of the proposed ASMC. 
- Case 2: the trifolium shape reference trajectory which is defined by (30): 


Xz = 5cos (= t) cos (= t) ‚Yz =5sin (= t) sin (= t) (30) 


The initial conditions are: [y,(0), yw (0), @(0)] = [5.1, — 1.2, 7| 


zl 


— reference trajectory f 
— -DOSM trajectory 
ASMC trajectory 


X(m) 


Figure 2. Trajectory tracking of the WMR in circular shape 


The simulation results are presented in Figures 4 and 5. Same as case 1, the trajectory of the WMR in 
both controllers are exact and robust against unknown input disturbances and the unknown wheel slips. But 
the tracking errors of ASMC controller converge to zero quickly in Figure 4 (about 0.5 s with ASMC 
meanwhile DOSM about 9 s). In Figure 5(a), at the time where the orientation of the WMR changes a lot, the 
trajectory controlled by the proposed controller has lower tracking error than the DOSM controller. In the 
x-axis, the maximum error tracking of the proposed controller is 0.1 m and the DOSM controller is about 
0.75 m (at the time the disturbances act on the system). Figure 5(b) is the comparison of the control torque 
between two control methods. It can be seen that the proposed control method is smoother, and the chattering 
problem is eliminated. Chattering reduction will increase the life of the actuator, this is an advantage of the 
proposed control method. Moreover, when the disturbances are acting on the system, the control torques of the 
DOSM exhibit an increase exceeding 40 Nm, whereas the proposed ASMC maintains these values below 20 Nm. 
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Figure 3. Comparative response of tracking error and motor torqure, (a) trajectory tracking error in X, Y 
direction and (b) comparison of the wheel motor torque 
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Figure 4. Trajectory tracking of the WMR in trifolium shape 
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Figure 5. Comparative response of tracking error and motor torqure, (a) trajectory tracking error in X, Y 
direction and (b) comparison of the wheel motor torque 


5. CONCLUSION 

The paper introduces an ASMC algorithm for a WMR considering wheel slips, uncertainties, and 
input disturbances. The state space model of the system has been constructed, and an adaptive SMC scheme is 
developed using this model to ensure system stability. Then an observer is designed to estimate the unknown 
components in the system which enhances the overall performance of the scheme. With such a design approach, 
this methodology can be effectively extended to a class of non-linear uncertain models as well. In this paper, 
the stability of the overall system was mathematically proved based on the Lyapunov theory. Simulation 
examples have been presented two types of reference trajectory. In each simulation, the responses of the 
proposed controller were compared with another SMC. The simulation results showed that chattering the 
performances of the proposed ASMC are significantly improved. 
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